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Abstract 



We study the Casimir piston for massless scalar fields obeying Dirichlet boundary conditions 
in a three dimensional cavity with sides of arbitrary lengths a. b and c where a is the plate 
separation. We obtain an exact expression for the Casimir force on the piston valid for 
any values of the three lengths. As in the electromagnetic case with perfect conductor 
conditions, we find that the Casimir force is negative (attractive) regardless of the values of 
a, b and c. Though cases exist where the interior contributes a positive (repulsive) Casimir 
force, the total Casimir force on the piston is negative when the exterior contribution is 
included. We also obtain an alternative expression for the Casimir force that is useful 
computationally when the plate separation a is large. 
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1 Introduction 



Two years ago, in an interesting paper P] , the Casimir piston was studied for a two-dimensional 
scalar field obeying Dirichlet boundary conditions on a rectangular region. Among other things, 
it was shown that the Casimir force on the piston is always attractive (negative) regardless of 
the ratio of the two sides. In this paper, we study the three-dimensional Casimir piston for 
massless scalar fields. A Casimir piston in three dimensions is depicted in Fig. 1. We choose 
the base to be a b x c rectangular region and a to be the plate separation (the distance from 
the base to the piston). The piston divides the volume into two regions. We refer to region I 
as the interior and region II as the exterior. Both regions contribute to the Casimir force on 
the piston. The Casimir piston therefore modifies some previous standard Casimir results [2] 
where the effects of the exterior region are not included. 

The Casimir piston for the electromagnetic field with perfect-conductor conditions in a three- 
dimensional rectangular cavity (box) was studied recently [3j and it was shown that the Casimir 
force on the piston is again attractive (in contrast to results without exterior region where the 
force could be positive). The piston for perfect-conductor conditions including the effects of 
temperature was studied further in [U [U [6] where among other things, the long and short 
distance behavior of the free energy was investigated. A theorem was obtained in [7], where 
it was shown that the Casimir force between two bodies related by reflection is always attrac- 
tive, independent of the exact form of the bodies or dielectric properties. This theorem was 
then generalized further in [8] where it was shown that reflection positivity implies that the 
force between any mirror pair of charge-conjugate probes of the quantum vacuum is attractive. 
Attraction does not occur in all Casimir piston scenarios. In a recent paper [9], the Casimir 
piston for a weakly reflecting dielectric was considered and it was shown that though attrac- 
tion occurred for small plate separation, this could switch to repulsion for sufficiently large 
separation. Moreover, for thick enough material, the force remained attractive for all plate 
separations in agreement with the results in [3] . Two recent preprints [lOj [11] also discuss 
scenarios where repulsive Casimir forces in pistons can be achieved. 

For the case of a massless scalar field in a three-dimensional cavity, approximate expressions 
for the Casimir force were obtained valid for small plate separation [3j. In this paper, we 
consider the general case of arbitrary lengths. We present exact expressions for the Casimir 
force on a piston due to a massless scalar field obeying Dirichlet boundary conditions in a 
three-dimensional box with sides of arbitrary lengths a, b and c. We find that the Casimir 
force on the piston is negative and runs from — oo (in the limit a — > 0) to (in the limit 
a — > oo). For small plate separation a, we recover the results found in [3]. We also obtain 
an exact alternative expression for the Casimir force that is useful computationally when the 
plate separation is large. We focus our attention on Dirichlet instead of Neumann boundary 
conditions because it is the more interesting case of the two. It is clear that Neumann boundary 
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conditions will yield a negative Casimir force since the contribution from both the interior and 
exterior are negative. It is not a priori obvious that in the Dirichlet case the Casimir force will 
be negative because there exists values of the ratios a/c and b/c where the interior contributes 
a positive (repulsive) Casimir force. It is therefore interesting to see that in such cases the 
exterior contributes a negative force of larger magnitude with the important consequence that 
the total Casimir force is negative. It is worth mentioning that the study of massless scalar 
fields is not only of theoretical interest but has direct relevance to physical systems such as 
Bose-Einstein condensates [121 E3 E] • 

The Casimir energy can be viewed as the energy with boundary conditions (a sum over discrete 
modes) minus the energy without boundary conditions (a volume integral over continuous 
modes). The sum over the discrete modes can typically be decomposed into a volume divergent 
term (the continuum part that can be subtracted), a surface divergent term and a finite 
part. In previous set-ups without region II, finite results were obtained by throwing out the 
surface divergent term. Though the finite results agreed with the zeta function regularization 
technique, there is nothing that can physically justify throwing out the surface term. It yields 
a cut-off dependent Casimir force that cannot be removed via a renormalization of the physical 
parameters of the theory |15^ [TBI [T7] . The agreement between zeta function regularization and 
cut-off technique (with surface term thrown out) occurs because the zeta function technique in 
effect renormalizes the surface term to zero. The Casimir piston resolves this issue satisfactorily 
by having the exterior and interior contributions to the surface divergence cancel. This has 
been demonstrated in Refs.[JJ [3] and we assume this cancellation to hold here. One can simply 
calculate the Casimir force F\ and F2 on the piston due to region I and II respectively without 
including the cut-off dependent terms. The total Casimir force on the piston can then be 
obtained by adding F\ and F%. One must just keep in mind that F\ and F2 actually have 
cut-off dependent terms but that they cancel when the two are added. 

There are two positive aspects to the Casimir piston: the exterior is now included in the 
calculation of the Casimir force (we add F2 to F\) and the surface divergence is handled via a 
cancellation procedure instead of simply throwing it out. 

We work in units where H = c = 1 (c is the speed of light). Note that from now on, when the 
variable c appears in the text, it always refers to one of the lengths of the base (see Fig. 1). 

2 Casimir piston in three dimensions: exact results 

The Casimir energy Ed for massless scalar fields in a d-dimensional box of arbitrary lengths 
L\, Ld obeying Dirichlet boundary conditions can be conveniently expressed as an analytical 
part - composed of Riemann zeta and gamma functions - plus a sum of over Bessel functions 
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(eq. QUID ; see appendix A and Refs.[T8| [23]): 



3=0 \ d) 

where Rj represents the sum over modified Bessel functions K u : 



^ = — m — ■ ( } 



n=l li=—oo 
i=l,...,j 



The prime in the sum for i?j means that the case where all ^'s are simultaneously zero is 
excluded. Note that Rj is a function of the ratios of the lengths. In (|2.ip . there is an implicit 
summation over the integers fej. The symbol C^" 1 is defined as 

= fl if k l <k 2 <...<h j ; l<kj <d-l 
ki,..,kj |q otherwise . 

The above symbol apparently does not have a name and we refer to it as the ordered symbol 
in Appendix A. The ordered symbol ensures that the implicit sum over the k{ in (|2.ip is over 
all distinct sets {k\, . . . ,kj}, where the ki are integers that can run from 1 to d — 1 inclusively 
under the constraint that k± <k 2 < ■ • • < kj. The superscript d— 1 specifies the maximum value 
of kj. For example, if j = 2 and d = 4 then ^j -1 fc . = Cj^ ^ 2 and the non-zero terms are £1^ 
, £1,3 and £2,3- This means the summation is over {k±,k 2 } = (1,2), (1,3) and (2,3). Note 
that the implicit summation over ki is also performed in Rj since Rj = Rj^L^/L^, .., Lj-./L^). 
For the special case of j = 0, Rj is defined to be zero and S,^ 1 k . and L k] are defined to be 

identically one so that ^f" 1 k . (l d y+i = f° r J = 0- 

From (|2.ip we can readily obtain the Dirichlet Casimir energy in three dimensions (d = 3): 
where i? is a function of L\,L 2 and L3 and represents the sums over Rj's i.e. 



7T 
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R(L 1 ,L 2 ,L 3 ) = T7rn[ L l R l {L l /L 3 ) + L 2 R 1 {L 2 /L 3 ) - R 2 (L 1 /L 3 , L 2 /L 3 ) (2.5) 



7T Li L 2 



where Ri{L\/L 3 ) means that R\ is a function of L\/L 3 . The functions i?i and R 2 are sums 
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over modified Bessel functions given by (|2.2p i.e. 



Ri(Li/L 3 



4n L 3 / Li 
^ L 3 



7T^ Li 

n=l 1=1 



R2{Li/Ls, L2/L3) 



tt n 

E E' 

n=l €i^2=— co 



(2.6) 
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L 3 / 



+ u 
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The Casimir energy does not depend on which sides are labeled L\, L2 and L3. Expression 
(|2.4j) for the Casimir energy is therefore invariant under permutations of the labels L\, L2 and 
L3 and we are free to label the three sides as we wish. For the Casimir piston depicted in 
Fig. 1, there are two regions to consider. In region I, the three sides are a, b and c and we 
label them L\ = c, L2 = b and L3 = o. In region II, the three sides are s — a, c and b and we 
label them L\ = s — a, L2 = c and L3 = b. The Dirichlet Casimir energy in region I and II is 
obtained by substituting the corresponding lengths in (|2.4p : 



Edi 
Ed2 



TT 2 be 



7T 



1440 a 3 + 32 vr a 2 
- 2 (s-a)c C(3) 



C(3) 



■(6 + c) 



7T 

96a 



+ R(c, b, a) 



TT 



(2.7) 



1440 6 3 

The function R(c,b,a) is obtained from (|2.5p and (|2.6p : 

7T c6 



327r62 ( S -« + c)- — + i?( S -a,c,6) 



i?(c, 0, a) 



TT 



16 a- 



■[ci?i(c/o) +bR 1 (b/c 



16 a 3 



i?2(c/a, fr/a) 



1 x x ' 

— J^Z^ J [^i(2vrn^c/a) + K^irnib/a) 



n=l 1=1 



be \ - 



00 n 3/2 x 3 ^ 2 ( 2vrn 



+ 



a J 



(2. 



71=1 



^c\ 2 + /^\ 2 
a / V a J 



1 3/4 



where the prime in the sum means that the case £1 = ^2 = is excluded from the sum The 
Casimir force on the piston is obtained by taking the derivative with respect to the plate 



1 Only the case when l\ and £2 are simultaneously zero is to be excluded. In particular, one can have £\ = 
when £2 7^ and vice versa. 
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separation a and then taking the limit s — > oo : 
d 



da 



[Edi + Em 



TT 2 bc C(3)(6 + C) 7T 7T 2 C C(3) 

' 480 a 4 + 16 7ra 3 96 a 2 ~ 1440 b 3 + 32 vr 6 2 



■R (c,b,a)— limR'(s — a,c,b) 



(2.9) 



We now evaluate the last two terms in (|2.9p . i?'(c, b,a) = R(c, b, a) can readily be obtained 
by taking the derivative of (|2.8p with respect to a: 

R'(c,b,a) = -7- YjYI 7 i K i fardel a) +K( (2vrn^/a) - -i^i (2vrn£c/o) - -.Ki fanlb/a 



n=l£=l 



(lis n 



n=l£ 1 ,£ 2 =— oo 



2 a (^ 2 c 2 +^ 2 6 2 



,3/4 



;£ 2 c 2 +£ 2 6 2 



,3/4 



} 



(2.10) 



where a prime on the Bessel functions denotes derivative with respect to the plate separation 
a. The last term in (|2.9|) can be written as 

d 

- lim — — R(u, c, b) 

u—>oo O U 



d 

lim R'(s — a, c, b) = lim — — R(s — a, c, b) 

s—foo s— >oo fj (X 



r 9 

m^oo a u 4 



^ oo oo 

"feJ^E 7 [^i( 27r ™^A0 +-Ki(2vrn^c/6) 



71=1 £=1 



863 



oo n 



3/2 



£ £' 

n=l £1/2=— oo 



(2.11) 



M'W(V)" + (x 



-1 3/4 



where the substitution u = s — a was made and R(u,c,b) was obtained from (12.80 by substi- 
tuting the appropriate lengths. The modified Bessel functions and their derivatives decrease 
exponentially fast so that the only term in (|2.1ip that survives is the case l\ = in the double 
sum. With i\ = 0, the remaining sum over £2 does not include zero and can be replaced by 
twice the sum from 1 to 00. One therefore obtains 



00 00 , z., I/O 

c \ — v /no\ i l I . 



limR'(s-a,c,b) = — 3 £ £ (-) K^ntc/b) 



71=1 £=1 



(2.12) 
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After substituting (|2.12p into (|2.9p . the Casimir force on the piston is 

ttHc C(3)(6 + c) TT 

b = 7 H 5 ^ — R (c, o, a) 

480 a 4 16 7ra 3 96 a 2 v ' ' ; 



vr 2 c C(3) c ~ ~ /„&N3tt ( 2 - 13 ) 

n=l £=1 



1440 b 3 32 it b 2 



Eq. (|2.13p is an exact expression for the Casimir force on the piston for Dirichlet boundary 
conditions. No approximations have been made. With R'(c,b,a) given by (|2.10p . one can 
calculate exactly the force for any values of a, b and c. Note that the second row in (I2.13P has 
no dependence on a and corresponds to the contribution from region II. If we set b = c and 
take the small a limit (a«6), we recover the expression for the Casimir force obtained in [3]. 
In this limit R'(c, b, a) is exponentially suppressed (exactly zero in the limit a — > 0) and with 
b = c, the second row in (I2.13P yields 0.004831546/c 2 in agreement with Dirichlet results in [3]. 

When a is sufficiently large, R'(c,b,a) dominates over the other a-dependent terms in (|2.13p . 
In fact, in the limit a — > oo, the other a-dependent terms vanish while R'(c,b,a) reduces to 
a finite function of b and c. Therefore, a full analysis of the Casimir force on the piston - 
one that goes beyond small values of a - requires one to have the exact expression (|2.10p for 
R'(c,b,a). 

In (|2.13p . the first and second row are the contributions from region I and II respectively: 

7T 2 bc C(3)(& + C) 7T 

= r H t; K — K (c, o, a) 

480 a 4 16 vr a 3 96 a 2 v ' ' ' 

TT 2 c C(3) c ^^/n6x3/2 ( 2 - 14 ) 

andF2 = "ll4W + ^-4l^£y *3 /2 (2Wc/ ). 

71=1 £=1 



To compute F\ and F2 we specify the two ratios a/c and b/c and express results in units of 
1/c 2 . Let us look at the case of the cube: a/c = 1 and b/c = 1. Using (|2.10p . we obtain 
R'(c,b,a) = —0.000214214. The last term in F2 - the sum over the Bessel function - yields 
—0.000271643. The remaining analytical terms in F\ and F2 can easily be evaluated. F\ and 
F2 for the case of the cube is given by 

F 1 . = -.005458275 +0.000214214 = -0.005244061 
F 2 . = 0.005103189 - 0.000271643 = 0.004831546 . 

We see that the Casimir force from region I is attractive and the force from region II is repulsive. 
Clearly, region II weakens significantly the total Casimir force. However, F2 is not large enough 
to reverse the sign and the Casimir force remains attractive: 

Fcube = F 1 + F 2 = -0.000412515 . (2.16) 
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The force Fi can actually be positive (repulsive) [SI] . For example, if a/c = 0.1 and b/c = 0.1 
then F\ = +3.80553076. However, the force due to the second region is then negative and larger 
in magnitude: F<i = —5.65818384. Adding the contribution from region II therefore causes a 
reversal of sign to take place. Though F\ is positive, the total Casimir force, F = F\ + F<i , is 
negative and equal to —1.85265308. 

The expression for the Casimir force on the piston, eq . (j2. 13[ ) . is valid for any positive values of 
a, b and c but is most useful computationally when the plate separation a is the smallest of the 
three lengths. The ratios b/a and c/a are then greater than or equal to one (we are also free to 
label the sides of the base such that c > b so that c/b\s also greater than or equal to one). The 
sums over the Bessel functions and their derivatives in (12,10p then converge exponentially fast 
yielding accurate and quick results. In Appendix B we derive an alternative expression F a n 
for the Casimir force on the piston that is useful computationally when the plate separation a 
is not the smallest of the three lengths. The alternative expression is given by (|B.7p : 



oo oo 



46 

n=l 1=1 



d_ 

da 



^ oo^ « - n 3 / 2 K 3/2 (2 irn^(^) + 

153 2_> 2_> 2_> 

n=l £i=l £2=— 00 



2 s (2.17) 



£io\2 /£ 2 c\ 2 



bJ + Kb 



3/4 



where the prime above the modified Bessel function K\ implies partial derivative with respect 
to a: K[(2Trn£a/b) = -^K\{2'K nla/b). As before, we are free to label the base such that 
c > b. If the plate separation a is not the smallest length, it follows that a >b and the above 
sums over Bessel functions and their derivatives converge exponentially fast. Both expressions, 
(|2.13p and (|2.17p . yield the same value for the Casimir force. However, computationally, 
expression (I2.13P is better to use if a is the smallest length and the alternative expression 
(|2.17p is better to use otherwise. 

For a given value of b and c, the Casimir force F on the piston ranges from —00 to zero 
corresponding to the two extreme limits of the plate separation a i.e. 

lim F = -00 and lim F = . (2.18) 

a^0 a— >oo 

The first limit in (I2.18|) follows readily if one uses expression (I2.13|) for the Casimir force. In 
the limit a— >0, the — 1/a 4 term dominates and goes to —00 (note that ]im a -+o R'(c,b, a) = 0). 
The second limit in (I2.18P follows readily if one uses the alternative expression F a u given by 
(|2,17p . In the limit a— >oo, F a i t is clearly zero since the Bessel functions and their derivatives 
decrease exponentially fast to zero as already mentioned at the end of Appendix B. One can 
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also understand this latter result intuitively: as a — > oo, region I becomes equivalent to region 
II and the forces from each region balance each other out i.e. lim^oo F\ = —F 2 . 

A plot of the Casimir force F versus a/c is shown in Fig. 2 for the case b/c = 1 (the force is 
in units of 1/c 2 ). The Casimir force is negative, has a large magnitude at small values of a/c 
and decreases rapidly in magnitude towards zero as a/c increases in agreement with the two 
limits given by (|2.18p . One obtains a similar plot for any value of b/c. A 3D plot of F versus 
a/c and b/c is shown in Fig. 3. The Casimir force is negative throughout and a slice taken at 
any value of b/c yields a similar profile to the 2D plot in Fig. 2 with the magnitude of the force 
shifting to greater values as b/c increases. For any given slice, the Casimir force lies between 
the two limits given by (12.181) . 



A Explicit expression for Dirichlet Casimir energy for <i-dimensional 
box with sides of arbitrary lengths 

In this appendix, we derive the explicit formula (|2.ip for the Casimir energy of massless scalar 
fields confined to a d-dimensional box of arbitrary lengths for Dirichlet boundary conditions. 
We begin by stating explicit formulas for the ci-dimensional Casimir energy obeying periodic 
boundary conditions [18] . The second step is to express the Dirichlet energy as a sum over 
the periodic energy The third step, the main part of the appendix, is to perform 

explicitly this sum to obtain the compact expression fQ2l) for the Dirichlet energy. 

The Casimir energy for massless scalar fields in a cf-dimensional box of arbitrary lengths 
Li, Ld and periodic boundary conditions can be explicitly expressed as an analytical part - 
composed of Riemann zeta and gamma functions - plus a sum of over Bessel functions [18] : 

V, M = " "E ft^k (r<¥) ^ CO + 2) + *,) 

j=o v i +1> (A.i) 
-7r C(3)£i ir 2 LiL 2 7rLi ttLiL 2 

— — n 5 r " * " — -ill 7T~ — -ilo 5 r * ' ' 

6L1 2vr L\ 90 L\ L\ h\ 

where Rj represents the sum over modified Bessel functions K v : 

~ ^,2^ K ^ ( 27T n a/^i-T^) 2 + ••• + & 7^) 2 ) 

R i = EEt — — m ■ (A ' 2) 

n=1 '=w [(^iSr) 2 + ■ ■ ■ + i^) 2 ] 4 

2 [E] uses a multidimensional cut-off technique and |23] uses the Epstein zeta function [19] technique. This 
technique has been developed extensively over the years |23]-[40| and there are some excellent books on the 
subject [201I2T1I22] . 
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The prime in the above sum means that the case where all €s are zero is excluded. Note that 
for j = one sets Rj to zero and Lj identically to one so that (j^.^'y+i is equal to 1/L% for 
j = 0. Note also that Rj is a function of ratios of lengths i.e. Rj = Rj(L\j Lj + i, Lj/Lj + i). 
The notation E p (d) is a compact way of saying that the Casimir energy E p is a function 

of the dimension d and the lengths Li, 

Our goal is to obtain a similar explicit expression for the case of Dirichlet boundary conditions. 
We begin by noting that the Dirichlet case can be expressed as a sum over the periodic Casimir 
energies E p (see [151125]): 

En = ^ £ (-D i+ ™ £ <«) 

m=l {ki,..,k m } 

ki<k 2 <-< 
k m <d 



The sum over the fcj's is over all sets {/ci, .., k m }, where the ki are integers that can run from 
1 to a maximum value of d under the constraint that k\ < &2 < ■ ■ ■ < k m . To specify that d is 
the maximum value we write k m < d under the sum in (|A.3|) . E p (m) is the periodic 

- L fcl'"' L '=m 

energy (TAJ]) repl acing d by m and L\ by L^, L2 by Lj~ 2 , etc. Note that the replacement L\ 
by Lfc i; etc. must also be performed inside Rj given by (|A.2|) . The above notation for the sum 
over ki is cumbersome. It is convenient to introduce a symbol £ % fcm defined by 

d fl if A;i<fc 2 <- ■ ■< k m ; 1 < ki < d 

I (J otherwise . 

The superscript d specifies the dimension which is the maximum value of k m . The above 
defined symbol apparently does not have a name and for simplicity we shall refer to it as the 
ordered symbol. Equation (|A.3h can now be conveniently expressed with the ordered symbol: 

E ° = WlJ2 ^r m ^L,k m E PLk ^ Lkm (m) (A.5) 

m=l 

where implicit summation over the k^s is assumed. After substituting (jA.ip into (|A.5[) one 
obtains 



771—1 



^ = E (-l)^t { E SP^t ( r (^) CO' + 2) + ^) } (A.6) 

m=l i=o v 

where i?j is the function (|A.2|) with Li replaced by L^, L2 by Lfc 2 , etc. For simplicity we 
define 

. >+ ^fe^( r(¥) ^ c<i+2)+ ^ (A ' 7) 
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and rewrite (|A.6P as 

d m—1 d—1 d 

£d = 2wE $Z( _1 ) d+m ^iv,fcm/i fel ,.., fcj+1 = ^TiE E ( _1 ) d+m+1 ^*i,»,*m^'fe ll ..^ +1 

m=l j=0 j=0 m=j+l 

(A.8) 

where we have rewritten the limits on each sum. We can decompose £j? km into a sum of 

two terms: Cf" 1 & d + ^fc^ 1 fcm' ^ e ^ rs ^ term, ^fc^" 1 fc <2 ' means that A; m is set to its 
maximum value of d and the sum is now over the remaining ki's with £; m _i having a maximum 
possible value of d — 1. With fc m = d in the first term, the maximum possible value of k m 
in the second term is d — 1 (hence the superscript — 1 in the second term). Note that for 
the special case m = d, the above decomposition yields only one term not two terms i.e. 
£fci,..,fcd = ZfaXkd-i,* + since can om y De equal to d. 
With this decomposition the sum over m becomes 

d d 

-d-i i 



E, i \rf+m+l _ \ / -i\d+m+l ft d—1 _i_ e d— 

l- X J ?fc X ,..,fe m - ^~ ^ L^fcl,..,fcm-l,d" , "^fcl, 

m=.j+l m=j'+l 

= f_1 W+J [ P d_1 -\-(fd-l _ td-l ) _ ( fd-l _pd-l 

\ / L ' k\,..,kj,d ~ \ ^ > / V s ki,..,kj+2 ^ ki,..,kj+2>o 



(AS 



The two terms inside each pair of round brackets in (1A.9|) have opposite signs and cancel each 
other H so that the sum over m reduces to only the first term 

E (-i) d+m+1 ,..,*„, = • (a.io) 

m=j+ 1 

The Dirichlet Casimir energy is obtained by substituting (lA.lOj) in (1A.8|) : 

The function fj k k . d ' ls obtained by setting kj+i equal to d in (IA.7jl . We finally obtain our 
explicit expression for the Dirichlet Casimir energy 

£o = E(-l)" +J «J-JU {%w 1 ( r (4 2 ) »^ CO + 2) + *,) } (A.12) 
j=o y dj 



'in the first pair of round brackets 1 k . 1 cancels with — 1 fe . l d . The fact that Ay+2 is equal to d in 



the latter term is irrelevant since the summation over fj in (|A.8|I stops at fc,-+i for a given j. Therefore 

1 1 • • • fc j . i 

iti 1 k +1 d ' s equivalent to ^J?" 1 fc. +1 - The same logic applies to the terms inside the other round brackets. 
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where Rj is given by (|A.2[) with L\ — » L^, — ► ^fc J+1 = -^d i.e. 



^ = — : ; — ~m — • (A - 13) 



n=l li=—oo 
i=l,...,j 
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For the case j = 0, Rj is zero and 1 fc . and are defined as unity. 



B Alternative expression for Casimir force on piston 

One can derive an alternative expression for the Casimir force F on the piston by labeling the 
lengths Li, L2 and L3 differently. We are free to label the lengths in any way we want since 
the Casimir energy is invariant under permutations of L±, L2 and L3. In region I, the three 
lengths are a, b and c and we label them now L\ = a, L2 = c and L3 = b. In region II, the 
three lengths are s — a, b and c and we label them now L\ = s — a, L2 = c and L3 = b. The 
Dirichlet Casimir energy in region I and II is then obtained via (|2.4p 

„ 7T 2 ac C(3) , . 7T 

E ™ = - 1440 V + 3^<° + C) " 966 + ^ C ' b) (R1) 

„ 7T 2 (s — a) C C(3) , 7T „, 

E ^ = -im-^ + ^ s - a + c) -Wb + R ^ s - a ^ b) 

where R(a, c, b) and -R(s — a, c, 6) are defined via (|2,5p and (|2.6p . The Casimir force F\ due to 
region I and F2 due to region II (with s — > 00) are 



(B.2) 



The total Casimir force .F = i 7 ! + F2 on the piston is then simply 

F = -R'(a, c, b) - lim R'(s - a, c, b) (B.3) 

where the prime denotes derivative with respect to the plate separation a. Note that the 
analytical terms - the Riemann zeta and gamma terms - have canceled. This is also what 
occurs in the two-dimensional Casimir piston (see [1]). The second term in (|B.3P has already 
been obtained and is given by (12.121) . The function R(a, c, b) can be obtained from the function 
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R(c, 6, a) given by (|2.8|) by replacing c with a, 6 with c and a with 6 i.e. 
1 x Xl 

R(a,c,b) = j^YjYj J [^i(2vrn£a/6) + K^ir ntc/b) 



n=l l=\ 



oo oo n 3 / 2 K 3 / 2 \ 2irn 



ac \ -> ^-v 

863 



ii a 



+ 



(B.4) 



n=l £i,£2=—oo 



I + I - 



3/4 



The derivative of R(a, c, b) with respect to the plate separation a is 

1 oo oo 

R'(a,c,b) = -J2J27 K '^ nia / b ) 



46 
9a 



n=l 1=1 



oo oo n 3 / 2 i^ 3/ / 2 ( 2 7rn 



8 63 2_> 2^ 



4a\2 /^ 2 cn2 



+ 



n=l l\,£i=— oo 



il a 



+ 



^2C\2 



3/4 



(B.5) 



With R'(a,c,b) given by (|B.5|) and lim^oo R'(s — a, c, 6) given by (|2.12p . the Casimir force 
JHH) yields 



1 oo oo 



71=1 £=1 



+ 



da 



a c 
86 s " 



x oo n 3 / 2 i^3/ 2 (2 7rn 



E E' 

n=l £ 1 ,£ 2 =—oo 



4a\2 /£ 2 CN2 



+ 



2 + /^ 2 C\2 



n3/4 



(B.6) 



c ~ ~ n 3 / 2 ^ 3/2 (2 7rn^ 2 c/6 
453 2^ L^i 



n=l £ 2 =1 



?2_C\3/2 

T 



The above expression has three terms and it can be simplified by noticing that the l\ = case 
in the second term cancels out with the last term. The Casimir force on the piston reduces to 
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the following final expression: 



-.0000 



n=l £=1 



+ 



d_ 

da 



ac 
463 



00 00 00 n 

EE E- 

n=l £ 1= ie 2 =-Qo 



3/^ 3/2 ( 2 .„J(^)% 



he 



fia\2 



13/4 



(B.7) 



The above is our alternative expression for the Casimir force on the piston. It is valid for 
any positive values of a, b and c but it is especially useful computationally when a is not 
the smallest of the three lengths. We are free to label the base such that c > 6. If a is not 
the smallest length, then the ratios a/6 and c/6 are both greater than or equal to one. This 
ensures that the sums over the Bessel functions and their derivatives in (1B.7|) will converge 
exponentially fast making computations easy and accurate. Note that the sum over l\ and 
the sum over i in (|B,7j) do not include zero. Therefore as a increases the Bessel functions and 
their derivatives will always decrease exponentially and reach zero in the limit a — > 00. The 
Casimir force on the piston is therefore zero in the limit a 



00. 
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Figure 1: Casimir piston in three dimensions. 
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Figure 2: Casimir force F versus a/c for the case b/c = 1 where a is the plate separation and 
b and c are the sides of the base. The force is in units of 1/c 2 . The force is large and negative 
at small values of a/c and remains negative with its magnitude decreasing quickly to zero as 
a/c increases. One obtains a similar plot for any value of b/c (see the 3D plot below) 
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Figure 3: 3D plot of Casimir force F versus a/c and b/c. The force is in units of 1/c 2 . For a 
given b/c, the profile is the same as in the 2D plot: the force is large and negative at small values 
of a/c and remains negative with its magnitude decreasing quickly to zero as a/c increases. 
The value of b/c shifts the magnitude of the force towards larger values as it increases. 
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